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Abstract 

The  systems  used  to  protect  eyes  and  sensors  from  frequency-agile  laser  weapons  must 
be  capable,  over  a  relatively  broad  range  of  frequencies,  of  absorbing,  refracting,  deflect¬ 
ing,  or  scattering  laser  radiation  of  high  intensity,  even  while  affording  high  transmission 
to  light  of  low  to  moderate  intensity.  The  design  of  such  devices  begins  with  an  accurate 
characterization  of  the  nonlinear  optical  response  of  the  candidate  active  materials.  In 
this  work,  the  relevant  material  parameters  can  be  extracted  from  experimentally  obtained 
data  only  by  indirect  means:  that  is,  the  actual  experimental  results  are  compared  to  the 
results  obtained  from  a  theoretical  model  employing  hypothetical  values  of  the  desired 
parameters,  and  the  parameters  of  the  model  are  varied  until  the  model  results  match  the 
experimental  data.  If  this  procedure  is  to  work,  the  model  must  capture  all  the  essential 
features  of  the  interaction  of  the  material  to  be  characterized  with  the  laser  beam  propa¬ 
gating  through  it.  We  review  the  assumptions  and  approximations  underlying  the 
commonly  used  models  for  nonlinear  absorption,  paying  particular  to  the  limits  of  appli¬ 
cability  that  these  impose. 
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1  Introduction 


Over  the  years,  the  nonlinear  optical  response  of  an  enormous  variety  of  materials  has  been 
studied,  the  absorptive  and  refractive  contributions  to  that  response  have  been  accurately  meas¬ 
ured  and  characterized,  and  the  underlying  physical  mechanisms  have  been  elucidated.  Most 
materials  exhibit  saturable  absorption:  that  is,  they  absorb  the  greatest  fraction  of  the  incident 
radiant  power  when  the  input  fluence  (total  delivered  energy  per  unit  area)  is  low,  with  the  ab¬ 
sorption  diminishing  as  the  fluence  increases.  Physically,  this  occurs  because  the  absorption 
cross-section  for  the  ground  state  is  greater  than  the  absorption  cross-section(s)  for  the  excited 
state(s).  The  observed  decrease  in  absorption  with  increasing  fluence  results  simply  from  the 
depletion  over  time  of  the  ground  state  population  as  an  ever-increasing  fraction  of  the  absorber 
molecules  is  excited. 

The  opposite  situation,  in  which  the  absorption  cross  section(s)  of  the  excited  state(s)  ex¬ 
ceeds  that  of  the  ground  state,  gives  rise  to  reverse  saturable  absorption  (RSA),  in  which  the 
absorption  increases  with  input  fluence.  This  enhancement  of  absorption  reflects  the  increase  in 
the  population  of  the  (more  strongly  absorbing)  excited  state(s)  relative  to  that  of  the  (less 
strongly  absorbing)  ground  state.  First  reported  in  the  dyes  sudanschwartz-B  and  indanthrone  by 
Guiliano  and  Hess  [1],  RSA  has  since  been  observed  in  a  variety  of  organic  and  organo-metallic 
molecules  and  complexes,  including  ion  tricobalt  clusters  [2],  King’s  complex  [3],  metallo- 
phthalocyanines  [4],  metallo-napthalocyanines  [4,  5],  fullerenes  (Cgo)  [6-9],  metallo- 
tetrabenzporphyrines  [10],  and  cubane-like  transition-metal  clusters  [11]. 

Both  types  of  absorptive  nonlinearity,  saturable  absorption  as  well  as  RSA,  offer  practical 
applications.  For  instance,  saturable  absorbers  are  commonly  used  as  passive  Q-switches  in 
pulsed  lasers.  Materials  exhibiting  RSA,  on  the  other  hand,  could  be  exploited  in  optical  limiting 
applications  such  as  eye  and  sensor  protection  systems,  which  require  high  transmission  of  low- 
intensity  radiation  but  high  absorption  at  high  input  [12,  13]. 

The  same  population  redistribution  responsible  for  absorptive  nonlinearities  also  gives  rise 
to  changes  in  the  refractive  index.  (One  might  expect  this,  since  the  leading-order  nonlinear  re¬ 
fractive  index  and  nonlinear  absorption  coefficient  are,  respectively,  the  real  and  imaginary  parts 
of  the  third-order  nonlinear  susceptibility  and  thus  are  linked  by  Kramers-Kronig  relations  [14].) 
A  full  characterization  of  the  nonlinear  optical  response  of  a  material  involves  the  measurement 
of  both  absorptive  and  refractive  nonlinearities.  These  two  quantities  are  typically  measured 
with  Z-scan  techniques  [15,  16].  In  a  Z-scan  experiment,  a  sample  of  material  is  moved  from  one 
side  of  focus  to  the  other  along  the  propagation  path  (Z)  of  a  tightly  focused  Gaussian  beam  and 
the  apertured  and  unapertured  transmittance  of  the  sample  is  measured  as  a  function  of  the  sam¬ 
ple’s  position  relative  to  focus.  From  the  apertured  transmittance  curve  (“closed  aperture  Z- 
scan”),  one  may  deduce  the  sign  and  magnitude  of  the  nonlinear  refractive  index;  one  then  uses 
the  unapertured  transmittance  data  (“open  aperture  Z-scan”)  to  infer  the  strength  of  the  nonlinear 
absorption. 

In  characterizing  optically  nonlinear  materials,  one  generally  extract  the  relevant  material 
parameters  from  experimentally  obtained  data  only  by  indirect  means:  that  is,  the  actual  experi¬ 
mental  results  are  compared  to  the  results  obtained  from  a  theoretical  model  employing 
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hypothetical  values  of  the  desired  parameters,  and  the  parameters  of  the  model  are  varied  until 
the  model  results  match  the  experimental  data.  If  this  procedure  is  to  work,  the  model  must  cap¬ 
ture  all  the  essential  features  of  the  interaction  of  the  material  to  be  characterized  with  the  laser 
beam  propagating  through  it. 

In  this  report,  we  review  the  existing  models  for  nonlinear  absorption,  paying  particular 
attention  to  the  assumptions  and  approximations  involved  in  each  model  and  to  the  limits  of  ap¬ 
plicability  that  these  impose.  We  begin  our  review  with  a  discussion  of  the  analytic  model 
presented  by  Wei  et  al.  in  one  of  the  earliest  papers  on  the  open  aperture  Z-scan  technique  [4]. 
Inconsistencies  in  the  assumptions  that  underlie  the  Wei  model  entrain  numerous  difficulties,  but 
these  problems  are  partially  offset  by  the  convenience  of  a  fully  analytic  solution.  Section  3 
treats  the  next  most  complicated  excited  state  absorption  model.  This  model,  though  still  very 
simple,  can  no  longer  be  fully  solved  by  purely  analytic  means  in  the  general  case.  However,  a 
complete  analytic  solution  is  possible  in  the  limit  of  low  fluence,  and  we  outline  this  solution  in 
Section  4.  This  is,  to  the  best  of  our  knowledge,  a  new  result  that  does  not  appear  in  the  litera¬ 
ture,  and  it  represents  a  definite  improvement  over  the  original  model  of  Wei  et  al.  In  Section  5, 
we  conclude  our  review  of  the  simplest  excited  state  absorption  models  with  a  comparative 
discussion  of  the  three  models  described  so  far.  In  Section  6,  we  turn  our  attention  to  the  so- 
called  five-band  model,  a  numerical  model  broadly  applicable  to  porphyrins  and  other 
conjugated  aromatic  compounds.  Section  7  introduces  a  commonly  seen  approximation  in  the 
five-band  model,  and  Sections  8  and  9  assess  the  impact  of  this  approximation.  Section  10 
describes  further  simplifications  of  the  five-band  model  that  are  applicable  in  particular  limits. 
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2  An  Analytic  Absorption  Model: 

Excited  State  Absorption  With  Constant  Ground  State  Absorption 


Among  the  earliest  theoretical  treatments  of  excited  state  absorption  is  the  following  sim¬ 
plified  model,  in  which  the  excited  state  is  populated  by  linear  absorption  from  the  ground  state 
[4].  In  this  model,  the  number  density  ri\(t)  of  molecules  in  the  excited  state  is  given  by 


dnx  _  a  j 
dt  hv 


(1) 


where  v  is  the  frequency  of  the  laser  radiation,  h  is  the  Planck  constant,  and  I  the  irradiance;  a , 
the  linear  absorption  coefficient,  is  a  constant.  The  irradiance  is  determined  by 

|^  =  -(a  +  CT1n1)/,  (2) 

oz 

where  a!  is  the  excited  state  absorption  cross-section.  Though  highly  problematic  from  a  physi¬ 
cal  standpoint,  the  assumption  of  a  purely  linear  ground  state  absorption  process,  completely 
characterized  by  the  constant  coefficient  a ,  does  possess  one  overwhelming  virtue:  the  resulting 
model,  given  by  Equations  (1)  and  (2),  can  be  solved  analytically.  One  proceeds  as  follows.  In¬ 
troducing  the  total  fluence, 


00 

F(r,z)=  ^I(r,z,t')dt' , 

—00 


(3) 


one  integrates  (1)  and  (2)  over  all  time  and  combines  the  results  to  obtain 


d±=-aF-^iF\ 
8z  2  hv 


(4) 


For  a  sufficiently  thin  sample  (one  whose  thickness  L  is  much  less  than  the  Rayleigh  range  of  the 
beam),  one  may  disregard  the  diffraction  .of  the  beam  as  it  propagates  through  the  sample  and  so 
perform  the  longitudinal  ( z )  integration  of  (4)  over  thickness  of  the  sample.  This  yields  a  relation 
between  the  total  fluence  entering  the  sample  at  a  radial  distance  r  from  the  center  of  the  beam 
and  the  exit  fluence  at  the  same  r.  In  the  case  in  which  the  input  fluence  results  from  a  principal¬ 
mode  Gaussian  beam,  this  relation  may  be  integrated  over  the  transverse  beam  profile  to  give  the 
total  transmitted  energy  or,  equivalently,  the  normalized  transmittance  T  /  TLinear  : 


ln[l  +  <7] 


T 

1  Linear 


(5) 


Here,  the  dimensionless  parameter  q  is  given  by 


acTj 

2  hv 


F0  Leff’ 


where  F0  is  the  total  fluence  on-axis,  and  Lejy  =  (1  -  exp[-aZ,])  /  a  is  the  effective  sample  thick¬ 
ness. 
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As  we  have  already  pointed  out,  the  most  glaring  inconsistency  in  this  simplified  model  is 
introduced  by  assumption  that  the  rate  of  growth,  al,  of  the  population  density  of  molecules  in 
the  excited  state  depends  on  time  only  indirectly,  through  the  time  dependence  of  the  irradiance. 
In  fact,  this  rate  must  also  depend  on  the  number  of  molecules  available  to  be  excited,  that  is,  on 
the  ground  state  population.  It  is  for  this  reason  that  the  analytic  model  cannot  account  for  satu¬ 
ration  effects  that  arise  from  depopulation  of  the  ground  state.  In  addition,  the  absence  from  (1) 
of  a  term  describing  the  nonradiative  decay  of  molecules  from  the  first  excited  state  back  to  the 
ground  state  carries  the  tacit  assumption  that  the  time  scale  for  this  decay  is  much  longer  than  the 
laser  pulse  width.  In  much  the  same  way,  the  absence  of  terms  for  the  radiative  upward  transi¬ 
tions  between  the  first  and  second  excited  states  and  for  the  nonradiative  decays  in  the  opposite 
direction  implies  a  tacit  assumption  that  the  lifetime  of  the  second  excited  state  is  negligibly 
short.  This  last  assumption  presents  no  problem,  as  a  rule.  On  the  other  hand,  neglecting  nonra¬ 
diative  depopulation  of  the  excited  state  is  valid  only  in  the  case  of  (1)  extremely  short  pulses  or 
(2)  an  effective  excited  state  of  triplet  spin  multiplicity.  Such  a  state  would  be  populated  by  a 
combined  process  consisting  of  a  radiative  transition  from  the  (singlet)  ground  state  to  a  singlet 
excited  state  followed  by  a  nonradiative  transition  to  the  effective  triplet  excited  state  on  a  time 
scale  that  is  shorter  than  or  comparable  to  the  pulse  width. 
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3  Incorporation  of  Ground  State  Saturation  Effects: 
Two-Band  Rate  Equation  Model 


The  seminal  excited  state  absorption  model  described  in  the  previous  section  falsely  as¬ 
sumes  that  the  rate  of  energy  absorption  by  active  molecules  in  the  ground  state  is  independent  of 
the  number  density  nG  of  such  molecules.  A  more  realistic  treatment  replaces  the  constant  a  by 

<jG  nG(t) ,  where  aG  is  the  ground  state  absorption  cross-section.  In  this  case,  (1)  and  (2)  be¬ 
come: 


dnx 
dt 


^  =  ^nGI 
^  hv  G 


~  =  -(CGWG  +  CTlnlV 

8z 


(6) 

(7) 


Like  (1),  Equation  (6)  neglects 

•  nonradiative  decays  from  the  excited  state  back  to  the  ground  state  (a  valid  approximation  so 
long  as  the  excited  state  lifetime  is  much  longer  than  the  pulse  width);  and 

•  radiative  transitions  from  the  first  excited  state  to  the  second  excited  state  and  nonradiative 
decays  from  the  second  excited  state  back  to  the  first  (a  valid  approximation,  provided  the 
downward  transition  is  so  fast  that  molecules  reaching  the  second  excited  level  decay  almost 
immediately  to  the  first  excited  state). 

Since  the  lifetime  of  second  excited  state  is  taken  to  be  so  short  that  this  state  is  always  essen¬ 
tially  empty,  the  populations  of  the  ground  state  and  first  excited  state  satisfy 

«g(0  +  «i(0  =  "o>  (8) 


where  n0  is  the  total  number  density  of  absorbing  molecules. 


The  ground  state  absorption  cross-section  aG  is  determined  by  extinction  measurements  at 
low  irradiance  levels  and  is  related  to  the  linear  absorption  coefficient  by  a  =  n0  aG .  The  quan¬ 
tity  Fs  =  hv/oG  defines  a  natural  scale  for  the  fluence;  we  expect  that  saturation  effects  will 
become  significant  at  fluence  levels  approaching  Fs.  For  this  reason,  one  frequently  refers  to  Fs 
as  the  “saturation  fluence.” 


Equations  (6),  (7),  and  (8)  maybe  integrated  with  respect  to  time  by  the  same  procedure 
employed  before.  The  result,  which  appears  in  [13],  is 

=  -«o  { aiF  +  (<*g  "  0  ~  exp[-F  / Fj) } ,  (9) 

OZ 

or,  expressed  in  terms  of  the  dimensionless  quantities  /  =  F  /  Fs,  x  =  a z,  and  <;  =  cjj  /  ctg  , 

This  longitudinal  integration,  in  contrast  to  the  one  in  (4),  cannot  be  performed  analytically. 
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4  Improved  Analytic  Model 

The  weak  saturation  limit  of  (9)  suggests  a  slight  improvement  in  the  seminal  excited  state 
absorption  model  described  in  Section  2.  At  low  fluence  levels  (  F  /  Fs  « 1 ),  one  can  expand 
the  exponential  function  on  the  right-hand  side  of  (9)  and  so  obtain 

si=.aF.^i^slF\  do) 

dz  2  hv 


where  we  neglect  terms  of  order  (A  /  F$  and  higher.  This  result  is  identical  to  (4),  except  that 
the  quantity  Gj  —  oq  in  (10)  plays  the  role  of  Cj  in  (4).  Thus,  the  result  (5)  for  the  normalized 
transmittance  must  hold  in  this  case,  with  the  parameter  q  replaced  by 


.»  -1 

= — rr- — 


2  hv 


-1 


F e 


-a  L. 


‘eff' 


To  the  best  of  our  knowledge,  this  result  is  new  and  does  not  appear  in  the  literature.  It  of¬ 
fers  reasonable  utility  for  fitting  Z-scan  data  obtained  with  picosecond  (or  faster)  pulses,  for 
which  fluence  levels  are  considerably  lower  than  for  nanosecond  pulses  and  for  which  the  two- 
band  model  is  generally  a  good  approximation.  This  is  discussed  in  more  detail  in  Section  10. 
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5  A  Comparison  of  the  Simplest  Excited  State  Absorption  Models 

For  reasons  that  we  presently  discuss,  one  expects  the  two-band  rate  equation  model  de¬ 
scribed  in  Section  3  to  fit  the  open  aperture  Z-scans  of  typical  organic  absorbers  reasonably  well, 
particularly  if  the  curves  in  question  are  obtained  using  a  picosecond  (or  faster)  laser.  For  pur¬ 
poses  of  comparison  with  the  two  analytic  models,  both  of  which  are  limited  by  even  more 
serious  assumptions,  we  take  the  predictions  of  the  two-band  model  to  be  “exact,”  as  indeed  they 
become  on  the  appropriate  time  scales. 

In  contrast,  the  analytic  “constant  a”  model  of  Wei  et  al.  [4],  reviewed  in  Section  2,  pre¬ 
sents  obvious  difficulties.  As  we  have  already  mentioned,  the  approximation  involved  in 
characterizing  ground  state  absorption  by  the  linear  absorption  coefficient  a,  a  constant,  com¬ 
pletely  precludes  saturable  absorption.  Every  Z-scan  curve  computed  with  this  model  will  have 
the  “valley”  shape  characteristic  of  RSA,  regardless  of  the  relative  values  of  the  ground  and  ex¬ 
cited  state  absorption  cross-sections.  This  is  because  the  constant  a  model,  in  effect,  assumes  that 
all  the  active  molecules  are  constantly  absorbing  from  the  ground  state,  no  matter  how  many  of 
these  are  simultaneously  ascribed  to  the  excited  state.  This  failure  to  account  for  the  depletion  of 
the  ground  state  population  as  the  excited  state  population  grows  also  leads  the  constant  a  model 
to  overstate  the  overall  absorption  of  an  RSA  material  with  a  given  value  of  the  excited  state 
cross-section  CTj  . 
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Figure  1.  Theoretical  Z-scan  Curves  Obtained  From  Various  Excited  State  Absorption  Models. 
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The  improved  model  described  in  Section  4  is  superior  to  the  constant  cl  model  in  both  respects. 
It  is  capable  of  modeling  saturable  absorption,  and  it  overstates  the  absorption  of  RSA  materials 
to  a  lesser  degree  than  the  constant  a  model.  Of  course,  the  superiority  of  the  improved  model  is 
expected  to  be  greatest  in  the  weak  saturation  limit. 

The  theoretical  Z-scan  curves  displayed  in  Figure  1  graphically  illustrate  these  points.  All 
curves  were  calculated  for  a  solution  containing  10  active  molecules  per  cubic  centimeter 
(equivalent  to  a  1.66-mM  concentration).  In  every  case,  the  ground  state  absorption  cross- 
section  aG  was  taken  to  be  10"17  cm2,  corresponding  to  a  linear  absorption  coefficient  a  of  10 
cm-1.  The  left-hand  pair  of  plots  shows  the  theoretical  curves  predicted  for  a  reverse  saturable 
absorber  having  an  excited  state  absorption  cross-section  crj  of  4  x  10“17  cm2,  while  the  two 
plots  on  the  right  show  the  curves  for  a  saturable  absorber  with  a  CTj  of  0.5  x  10~17  cm2.  As  ex¬ 
pected,  the  constant  a  model  (long  dashes)  produces  the  same  valley-shaped  Z-scan  character¬ 
istic  of  RSA,  regardless  of  whether  the  material  is  saturably  absorbing  or  reverse  saturably  ab¬ 
sorbing. 

The  upper  two  plots  in  Figure  1  depict  low-fluence  Z-scans  such  as  might  be  obtained  from 
a  532-nm,  principal-mode  Gaussian  beam  delivering  100  nJ/pulse.  Such  a  beam  provides  an  on- 
axis  fluence  of  about  27  mJ/cm2.  Since  this  value  never  exceeds  74%  of  the  saturation  fluence  of 
37.3  mJ/cm2,  we  would  expect  the  approximation  (F / Fs  « 1)  that  gave  rise  to  the  improved 

model  to  hold  fairly  well.  In  this  case,  the  curves  generated  by  the  improved  model  (short 
dashes)  do  indeed  show  good  agreement  with  the  “exact”  curves  obtained  from  the  two-band 
model  (solid  lines)  for  both  saturably  absorbing  (top  right  plot)  and  reverse  saturably  absorbing 
materials  (top  left  plot).  On  the  other  hand,  the  lower  two  plots  illustrate  the  various  models’ 
predictions  for  Z-scans  at  fluence  levels  corresponding  to  an  energy  of  2  pJ/pulse,  which  is  fairly 
typically  for  a  Z-scan  performed  with  a  nanosecond  laser.  Assuming  a  spatial  beam  shape  that  is 
principal-mode  Gaussian,  such  a  pulse  provides  a  fluence  on  axis  of  about  550  mJ/cm  .  This  ex¬ 
ceeds  the  saturation  fluence  by  a  factor  of  almost  15,  so  our  expectations  for  the  improved  model 
should  not  be  particularly  high. 
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6  The  Five-Band  Model 


Nonlinear  absorption  in  a  wide  variety  of  organic  materials,  including  the  metallo- 
phthalocyanines  and  metallo-naphthalocyanines  which  have  been  the  focus  of  so  much  recent 
attention,  is  well  described  by  the  five-band  model  depicted  schematically  in  Figure  2.  Each 
molecule  is  assumed  to  lie  in  the  vibration-rotation  manifold  of  one  of  five  electronic  states:  the 
ground  state,  So;  one  of  two  excited  singlet  states,  S\  or  S2;  or  one  of  two  excited  triplet  states,  T\ 
or  72.  Let  ctg,  as,  and  a 7-  denote,  respectively,  the  absorption  cross  sections  for  the  transitions  So 
—>  S\,  S\  — >  1S2,  and  T\  — >  T2,  and  let  ks,  ks2,  kj,  kj2,  and  kisc  be  the  rate  constants  for  the  decays  S\ 
— >  So,  S2  — ^  Si,  T\  — >  Sq,  T2  — >  T j ,  and  S\  — >  1\ . 


Figure  2.  Five-band  Model. 


The  following  rate  equations  specify  the  time  evolution  of  the  number  densities  of  molecules  in 
the  So,  Si,  S2,  T\,  and  T2  bands;  we  refer  to  these  number  densities  as  “populations”  and  denote 
them  by  no,  ns,  ns2,  nr,  and  n-n,  respectively. 

~ — 7~nQ  I  +  ksns  +  kTnT  (11) 

ot  hv 

=  ^r~nG  1  -  (ks  +  kisc)nS  -  7^%  1  +  kS2”S2  (!2) 

ot  hv  hv 
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3hy 

dt 


hcns 


—  kjVij  — 


(J  rj-l 

—  nj’  I  +  kj2nT2 
hv 


(13) 


^  =  ^-„sl-ks2ns2  (14) 

dt  hv 

^Il=  £ LnTl_kT2riT2  (15) 

dt  hv 

As  before,  v  is  the  frequency  of  the  laser  radiation,  h  is  the  Planck  constant,  and  I  the  irradiance. 
Summing  these  equations,  one  verifies  that  they  are  consistent  with  the  condition  that  no,  the  to¬ 
tal  number  of  molecules  per  unit  volume, 

w0  =  nG  +  nS  +  nT  +  nS2  +  nT2  > 

be  constant  in  time.  Equation  (16)  may  be  used  to  eliminate  one  of  the  number  densities  (nn, 
say)  from  Equations  (11)  through  (15).  The  rate  equations  for  the  five-band  model  thus  consist 
of  a  system  of  four  coupled  ordinary  differential  equations  in  four  unknowns. 


We  stress  the  fact  that  Equations  (11)  through  (15)  describe  a  “five-band  model”  as  op¬ 
posed  to  a  “five-level  model.”  This  is  because  radiative  transitions,  such  as  the  singlet-singlet 
transition  from  the  ground  state  to  the  S\  band,  usually  result  in  excitation  high  into  the  vibration- 
rotation  manifold  of  the  target  electronic  band,  followed  by  a  rapid  decay  to  a  state  at  the  bottom 
of  the  band.  This  latter  state  is  nonresonant  with  the  incident  radiation.  As  a  result,  stimulated 
emission  terms  are  absent  from  (11)  through  (15). 


Transitions  between  singlet  and  triplet  states  are  forbidden  by  spin  selection  rules;  the  life¬ 
time  xT  =  kT  of  the  lowest  lying  triplet  state  T\  is  thus  extremely  long,  on  the  order  of 
microseconds.  In  contrast,  the  effects  of  spin-orbit  coupling  induced  by  the  metal  substitute  in 
metallo-phthalocyanines  and  metallo-naphthalocyanines  combine  with  the  closeness  in  energy  of 
the  S\  and  T\  bands  to  shorten  the  decay  time  xisc  =  kisc~}  for  the  intersystem  crossing  to  several 
nanoseconds. 


In  certain  organic  compounds  such  as  fullerenes  [8,  9,  17]  and  cubanelike  transition-metal 
clusters  [11],  there  is  evidence  that  the  dominant  mechanism  by  which  the  triplet  manifold  is 
populated  is  not  intersystem  crossing  but  rather,  a  process  of  ionization  and  subsequent  germi¬ 
nate  recombination.  In  this  case,  the  band  labeled  S2  in  Figure  2  actually  represents  an  ionized 
state  of  the  molecule.  The  following  modifications  are  required  to  adapt  the  rate  Equations  (1 1) 
through  (15)  to  systems  in  which  the  triplet  manifold  is  populated  primarily  by  ionization  fol¬ 
lowed  by  germinate  recombination: 

•  One  replaces  the  last  term  on  the  right-hand  side  of  (12),  kS2nS2 ,  by  ns  kion_e  nS2  . 

•  One  replaces  the  last  term  on  the  right-hand  side  of  (14),  ~kS2nS2  by 

—  (j^S  ^ T  )^" ion-e  ^ S2  • 

•  An  additional  term  nT  kion_e  n$ 2  is  added  to  the  right-hand  side  of  (13). 

Here,  kion  e  is  the  rate  constant  for  cation-electron  recombination  and  ns  and  nT  are  the  probabili¬ 
ties  of  re-forming  the  molecule  in  the  -Si  and  T\  state,  respectively.  (The  probability  of  re- 
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forming  molecules  in  the  ground  state  by  germinate  recombination  is  so  small  that  it  may  be 
safely  ignored  [17].) 

Completing  the  five-band  model  is  the  following  extinction  law,  which  describes  the  de¬ 
crease  in  intensity  of  the  propagating  beam  as  a  result  of  linear  and  nonlinear  absorption: 

—  =  —  +  ^SnS  +  •  07) 

oz 
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7  The  Truncated  Five-Band  Model 


One  can  frequently  simplify  the  five-band  model  by  assuming  that  the  populations  of  the  S2 
and  T2  bands  remain  at  a  constant  and  negligible  level  [13].  Except  at  extremely  high  irradian- 
ces,  this  approximation  should  be  generally  valid.  It  represents  an  enormous  simplification, 
allowing  us  to  neglect  the  last  two  terms  on  the  right-hand  sides  of  (12),  (13),  and  (16)  and  to 
ignore  Equations  (14)  and  (15)  entirely.  The  resulting  rate  equations  are 

=  -^-nG  I  +  ksns  +  krnT  (18) 

dt  hv 

^~  =  ^~nGI  ~(ks +kisc)ns  (19) 

at  hv 

dJhL  =  klscns-kTnT  (20) 

ct 

Using  the  following  molecule  conservation  equation, 

n0  =nG+nS+nT’  @1) 

one  may  eliminate  one  of  the  number  densities  from  Equations  (18)  through  (20).  The  irradiance 
is  described  by  (17),  the  same  extinction  law  as  the  five-band  model.  We  refer  to  the  model  de¬ 
fined  by  (17)  through  (21)  as  the  truncated five -band  model. 
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8  Validity  of  Second  Excited  Band  Truncation 

In  order  to  examine  the  limitations  imposed  by  the  approximation  described  in  the  preced¬ 
ing  section,  we  use  the  full  set  of  rate  equations,  (11)  through  (15),  to  compute  the  evolution  in 
time  of  the  relative  number  densities  of  molecules  in  each  of  the  five  bands.  Figure  3  shows  the 
number  densities  on  the  input  face  of  a  sample  of  silicon  naphthalocyanine  (SiNc)  illuminated  by 
a  laser  pulse  of  Gaussian  temporal  profile  and  full  width  at  half  maximum  (FWHM)  pulse  width 
100  ps  (upper  pair  of  graphs)  and  10  ns  (lower  graph).  We  performed  these  computations  by  tak¬ 
ing  ksi  =  kn  =  (1  ps)-1,  kf  =  (300  ns)-1,  and  employing  newly  published  values  for  the  remaining 
rate  constants  and  the  absorption  cross-sections  [18].  At  these  peak  intensities  (2  GW/cm2  for 
the  100-ps  pulse  and  1  GW/cm2  for  the  10-ns  pulse,  corresponding  respectively  to  fluences  of  0.2 
J/cm2  and  10  J/cm2  on  axis),  the  population  of  the  T2  band  is  not  negligible  relative  to  that  of  the 
T\  band,  nor  is  the  S2  population  negligible  relative  to  the  S\  population. 


10  20  30  40 

Time  (ns) 


Figure  3.  Band  Population  Fractions  at  the  Input  Face  of  a  Sample  of  SiNc. 


9  Full  and  Truncated  Five-Band  Models  Compared: 

Excited  State  Saturation 

Because  it  takes  the  lifetimes  of  the  second  excited  states  to  be  vanishingly  small  and  thus 
precludes  excited  state  saturation,  the  truncated  model  overstates  excited  state  absorption.  In 
principle,  the  values  of  excited  state  absorption  cross-sections  extracted  from  open  aperture  Z- 
scan  data  via  the  truncated  model  will  lie  below  their  true  values.  In  practice,  this  discrepancy  is 
completely  negligible,  since  Z-scan  experiments  are  typically  performed  at  such  low  intensities 
that  the  normalized  transmittance  computed  with  the  truncated  model  differs  only  slightly  from 
that  calculated  with  the  exact  model,  the  difference  being  most  pronounced  for  the  portion  of  the 
curve  near  the  focal  plane.  This  situation  is  well  illustrated  in  Figure  4.  Shown  are  the  theoreti¬ 
cal  Z-scan  curves  for  a  1.89-mM  SiNc  sample,  0.1  mm  thick.  The  upper  curve  is  calculated  with 
the  full  five-band  model  and  the  lower  curve,  with  the  truncated  model.  Both  computations  as¬ 
sume  a  laser  source  at  532  nm  delivering  7.25  pJ  per  7.12-ns  FWHM  Gaussian  pulse.  The  beam 
is  assumed  to  be  a  principal-mode  (TEMo,o)  Gaussian  beam,  focused  to  a  spot  of  radius  10.8  pm. 


Figure  4.  Z-scan  Curves  Obtained  with  Full  and  Truncated  Five-band  Models. 


In  modeling  power  limiting,  however,  the  situation  is  very  different,  as  Figure  5  illustrates. 
The  power  limiting  curves  shown  in  Figure  5  are  for  a  0.1 -mm  sample  of  SiNc.  The  upper  solid 
curve  was  calculated  with  the  full  five-band  model,  while  the  lower  solid  curve  employed  the 
truncated  five-band  model.  The  dotted  curve  shown  in  the  plot  on  the  right  extrapolates  low- 
energy  transmission.  In  this  case,  the  “exact”  full  five-band  model  predicts  a  decline  in  excited 
state  absorption  resulting  from  excited  state  saturation.  Of  course,  one  might  argue  that  at  the 
fluence  levels  where  the  discrepancy  between  the  two  models  becomes  significant,  modeling  ef¬ 
forts  are  already  problematic  due  to  the  onset  of  a  host  of  processes  that  are  not  yet  well 
quantified,  such  as  plasma  and  bubble  formation. 
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Figure  5.  Power  Limiting  Curves  for  the  Full  and  Truncated  Five-band  Models. 


The  computational  savings  that  may  be  realized  by  employing  the  truncated  version  of  the 
five-band  model  instead  of  the  full  version  are  less  significant  than  one  might  think.  Displayed 
in  Table  1  are  some  representative  results  for  the  normalized  transmittance,  obtained  during  the 
calculation  of  Z-scan  curves  just  discussed,  along  with  the  computation  time  required  for  each. 


Table  1.  Relative  Computation  Times  for  the  Full  and  Truncated  Five-band  Models 


Full 

Truncated 

z 

T 

Time 

T 

Time 

Zr 

Tl 

00 

Tl 

00 

0 

0.54159 

3177 

0.52661 

3238 

1 

0.59610 

3559 

0.58772 

2935 

5 

0.87425 

3092 

0.87086 

2296 

10 

0.95993 

3072 

0.95645 

2227 
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10  Two-  and  Three-Band  Effective  Models 

For  the  appropriate  time  scales,  the  five-band  rate  equation  model  can  be  simplified  by  the 
omission  of  an  entire  manifold  of  electronic  states.  For  instance,  in  the  case  of  picosecond 
pulses,  the  triplet  states  are  always  essentially  empty  and  the  system  is  well  described  by  a  three- 
band  model  with  a  single  excited  state  with  absorption  cross  section  a,  =  as.  Conversely,  for 
long  pulses  (pulse  widths  x pulse  »  zisc )  and  continuous  wave  beams,  one  may  choose  to  ignore 
the  effects  of  the  S,  ->  S2  transition,  particularly  when  aT  »  <75,  as  is  frequently  the  case.  In 
such  instances,  one  often  adopts  a  three-band  model  with  cxx  =  $oT,  where  the  triplet  yield  <|>  is 
the  fraction  of  the  excited  singlet  population  which  reaches  the  triplet  state.  And  even  when  both 
the  singlet-singlet  and  triplet-triplet  excited  state  transitions  are  of  comparable  importance,  it  is 
always  possible  to  define  an  effective  excited  state  cross  section  in  a  quasi  three-band  model 
valid  for  a  particular  pulse  width  and  fluence  [13].  We  wish  to  consider  a  generic  situation  in 
which  there  is  a  ground  state  and  one,  two,  or  more  excited  bands  of  unspecified  spin  multiplic¬ 
ity.  Figure  6  depicts  the  situation  with  a  pair  of  excited  bands. 


Figure  6.  Three-band  Effective  Model. 


The  arguments  from  our  discussion  of  the  truncated  five-band  model  in  Section  7  apply 
equally  well  here.  If  the  lifetime  of  the  second  excited  band  is  sufficiently  short  (the  usual  situa¬ 
tion  unless  the  first  and  second  excited  bands  differ  in  spin  multiplicity),  it  is  generally  necessary 
to  track  the  population  of  this  band  only  in  situations  involving  extremely  high  fluence  levels.  In 
all  other  cases,  we  are  well  served  by  the  following  truncated  system  of  rate  equations: 
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dnG  aG  r  r 

-z r  =  — rrnGI  +  ki  "i 

ot  hv 

dnes  _<*G„  j  jl  „ 

OT  «V 

where  &i  is  rate  constant  for  decay  from  the  excited  state  to  the  ground  state.  These  two  equa¬ 
tions  are  equivalent,  as  can  be  easily  seen  from  the  condition  (21)  for  conservation  of  absorber 
molecules,  which  in  this  case  reads  simply  n0  =  nG  +ni.  The  level  populations  in  the  resulting 

two-band  effective  model  are  thus  described  by  a  single  differential  equation,  which  is  conven¬ 
iently  expressed  in  terms  of  the  fraction  X  of  molecules  in  the  excited  state  ( X=n\  /  no): 

F)Y  J 

^-  =  -~0-X)-klX.  (22) 

St  F, 

Here,  Fs  is  the  saturation  fluence,  defined  in  Section  3  above,  as  Fs  =  hv/<jG ■  The  differential 
equation  for  the  irradiance  is 

^  =  -»0  Ml -*)  +  <*!*}/•  (23) 

cz 

The  system  of  coupled  equations  (22)  and  (23)  constituting  the  two-band  model  may  be  solved 
analytically  in  two  special  cases,  which  we  now  review  in  turn. 


10.1  Negligible  Excited  State  Decay 

If  the  excited  state  lifetime  ii  is  much  greater  than  the  pulse  width  (as  is  the  case  for  nano¬ 
second  pulses  when  the  excited  state  is  a  triplet  state),  then  ii  may  be  taken  to  be  effectively 
infinite  and  k\  =  n  -1  =  0.  In  this  case,  the  two-band  model  reduces  to  the  model  developed  in 
Section  3.  Equation  (22)  may  be  integrated  with  respect  to  time  to  give 

jr(0  =  l-exp[-F(/)/FJ] 


where 


t 

F{t)  =  \l{t')dt' 

—00 

is  the  fluence  at  time  t.  Substituting  the  above  expression  for  X(t)  in  (23)  and  integrating  the  re¬ 
sult  over  all  time  yields  (9). 


10.2  No  Saturation 

The  other  special  case  in  which  (22)  and  (23)  may  be  solved  analytically,  while  arising  far 
less  frequently  in  practice  than  the  case  of  negligible  excited  state  decay,  is  nonetheless  applica¬ 
ble  to  a  number  of  materials,  notably  indoaniline.  It  occurs  when  the  ground  and  excited  state 
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absorption  coefficients  are  identical:  gg  =  =  g.  In  this  case,  the  overall  absorption  coefficient 

is  strictly  constant, 

a  =  ®GnG  +alwl  =  a(nG  +nl)  =  an0’  (24) 

so  the  irradiance  is  independent  of  time  and  described  by  the  familiar  Beer’s  Law  relation: 

7(z)  =  /(z0)e-a(z-z«). 


Thus,  (22)  may  be  integrated  to  give 

1  -  exp 

m=— 


f,  /(*)! 

1  -  exp 

fi  ,  /(z)l 

t 

k,  + 

1  F 

V  rs  ) 

t 

1  + - 

V  J 

Ti. 

1  + 


/(z) 


1  +  -7' 


/(z) 


where  ti  =  k\  1  is  the  excited  state  lifetime  and 


h  =  *1  Fs 


hv 


GTj 


The  constant  Is  is  usually  referred  to  as  the  “saturation  intensity.”  (This  designation  is  something 
of  a  misnomer  in  the  present  case,  since  the  overall  absorption  exhibits  no  saturation!) 
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11  Conclusions 


This  report  summarizes  the  most  commonly  seen  models  for  excited  state  absorption.  In 
these  models,  nonlinearity  arises  from  the  redistribution  of  absorber  molecules  among  various 
available  energy  bands,  the  absorption  cross-sections  of  which  differ  in  magnitude.  The  existing 
analytic  model  for  excited  state  absorption  involves  approximations  that  entirely  preclude  satur¬ 
able  absorption  and  lead  one  to  underestimate  the  absorption  cross-sections  of  reverse  saturable 
absorbers  based  on  their  Z-scans.  We  present  here  an  improved  analytic  model  that  is  able  to 
reproduce  the  behavior  of  saturable  absorbers  and  that  mitigates  the  degree  to  which  model  fits 
of  Z-scan  curves  underestimate  the  excited  state  absorption  cross-sections  of  reverse  saturable 
absorbers. 

The  most  complete  description  of  the  conjugated  organic  compounds  currently  under  study 
as  potential  optical  limiters  is  provided  by  a  model  involving  five  molecular  energy  bands  and 
which  we  straightforwardly  designate  the  full  five-band  model.  One  frequently  assumes  that  the 
populations  of  the  highest  energy  band  of  each  spin  multiplicity  (i.e.,  the  S2  and  T2  bands)  remain 
at  a  constant  and  negligible  level.  This  approximation  gives  rise  to  a  model  that  we  refer  to  as 
the  truncated  five-band  model.  Most  of  the  theoretical  calculations  reported  in  the  recent  litera¬ 
ture  employ  this  latter  model,  which  most  authors  refer  to  simply  as  the  “five-band  model.”  We 
demonstrate  in  this  report  that,  at  least  at  the  center  of  the  beam,  the  populations  of  the  ignored 
bands  are  in  fact  non-negligible  even  at  moderate  fluence  levels.  These  populations  should  defi¬ 
nitely  be  included  (through  the  use  of  the  full  five-band  model)  in  the  modeling  optical  limiter 
performance,  though  they  may  be  safely  ignored  (by  the  use  of  the  truncated  model)  in  the  inter¬ 
pretation  of  Z-scan  data. 

In  modeling  the  propagation  of  nanosecond  pulses,  one  is  generally  restricted  to  the  use  of 
either  the  full  or  truncated  five-band  model,  assuming  that  the  material  to  be  modeled  has  an  en¬ 
ergy  level  structure  consistent  with  the  five-band  scheme.  For  very  short  or  very  long  pulses, 
additional  simplifications  are  possible,  and  these  give  rise  to  two-  and  three-band  effective  mod¬ 
els.  An  effective  model  may  also  prove  useful  in  a  preliminary  attempt  to  elucidate  the  nonlinear 
optical  behavior  of  materials  having  energy  level  structures  that  differ  from  the  standard  five- 
band  arrangement. 
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